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We study several properties of distillation protocols to purify multilevel qubit states (qudits) 
when applied to a certain family of initial mixed bipartite states. We find that it is possible to use 
qudits states to increase the stability region obtained with the fiow equations to distill qubits. In 
particular, for qutrits we get the phase diagram of the distillation process with a rich structure of 
fixed points. We investigate the large-D limit of qudits protocols and find an analytical solution in 
the continuum limit. The general solution of the distillation recursion relations is presented in an 
appendix. We stress the notion of weight amplification for distillation protocols as opposed to the 
quantum amplitude amplification that appears in the Grover algorithm. Likewise, we investigate 
the relations between quantum distillation and quantum renormalization processes. 



PACS numbers: 03.67.-a, 03.67.Lx 



I. INTRODUCTION 

The experimental analysis of the intriguing proper- 
ties of entanglement in quantum mechanics requires the 
availability of stable sources of entanglement. Despite 
the nice properties exhibited by entanglement, it has the 
odd behaviour of degrading by the unavoidable contact 
with the external environment. Thus, for the entangle- 
ment to be assessed as a precious mean, we must devise 
some method to pump it up to the entanglement source 
in order to sustain a prescribed degree of entanglement 
that we may need whether for quantum communication 
protocols (teleportation, cryptography, dense coding) or 
quantum computing (algorithmics) (for a review see Q, 
and references therein). 

Quantum distillation or purification protocols are pre- 
cisely those methods, that have been devised to regen- 
erate entanglement leakages of an entanglement source. 
Here we are interested in the purification of mixed states 
of bipartite type, having in mind the realization of a com- 
munication protocol by two parties, Alice and Bob. The 
seminal work of Q has provided us with a standard dis- 
tillation method that has been the focus to developing 
more protocols with the aim at improving its original 
performances. We shall refer to this distillation protocol 
as the BBPSSW protocol. There are feasible experimen- 
tal proposals for this type of protocols using polarization 
beam splitters (PBS) Likewise, there also exist meth- 
ods for the distillation of pure states that have been 
implemented experimentally 

In addition to the initial purpose for which the quan- 
tum distillation protocols were devised, they have found 
another very important application in connection to the 
problem of quantum error correction: quantum informa- 
tion needs to be protected from errors even more than 
classical information due to its tendency to become de- 
coherent. To avoid these errors, one can resort to the 
ideas of quantum error correction codes , [1| and fault- 
tolerant quantum computation 0, j^. However, entan- 
glement purification is another alternative to decoherence 



which gives a more powerful way of dealing with errors 
in quantum communication 

In a typical quantum communication experiment, Alice 
and Bob are two spatially separated parties sharing pairs 
of entangled qubits. The type of operations allowed on 
these qubits arc denoted as LOCC (Local Operations and 
Classical Communication): they comprise local unitary 
operators Ua 'Si Ub on each side, local quantum measure- 
ments and communication of the measurement results 
through a classical channel. These local quantum op- 
erations will suffer from imperfections producing local 
errors. Futhermore, Alice and Bob will also face trans- 
mission errors in their quantum channels due to dissipa- 
tion and noise. To overcome these difficulties, they will 
have to set up an entanglement purification method. In 
short, a protocol like the BBPSSW creates a reduced set 
of maximally entangled pairs (within a certain accuracy) 
out of a larger set of imperfectly entangled pairs: entan- 
glement is created at the expense of wasting extra pairs. 
The degree of purity of a mixed entangled pair is mea- 
sured in terms of its fidelity with respect to a maximally 
entangled pure pair, which is the focus of the purifica- 
tion protocol. After the BBPSSW protocol, a new dis- 
tillation protocol was introduced in |l2l | by the name of 
Quantum Privacy Amplification (QPA) which converges 
much faster to the desired fidelity [l3l. UM. Other pro- 



tocols known as quantum repeaters |l5j. [Til allow us 
to stablish quantum communication over long distances 
by avoiding absortion or depolarization errors that scale 
exponentially with the length of the quantum channel. 

The advantages of dealing with £)-dimensional or mul- 
tilevel quantum states (qudits) instead of qubits are quite 
apparent: an increase in the information flux through 
the communication channels that could speed up quan- 
tum cryptography, etc. [23, [IJ, Thus, it has 
been quite natural to propose extensions of the purifi- 
cation protocols for qudits. One of the proposals 0| 
relies on an extension of the CNOT gate that is unitary, 
but not Hermitian. Recently, another very nice proposal 
has been introduced ^3 based on a generalization 
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of the CNOT gate that is both unitary and Hcrmitian 
and gives a higher convergence. In this paper we make a 
study of the new purification protocols of [T^ , when 
they are appHed to mixed bipartite states of qudits that 
are not of the Werner form. In this way, we combine 
some of the tools employed by the QPA protocols |0| 
with the advantages of the new methods. 

This paper is organized as follows: in Sect, ^we re- 
view simple distillation protocols for qubits not in Werner 
states and we generalize them for the purification of any 
of the Bell states. In Scct. lIIII wc extend the previous pro- 
tocols to deal with multilevel qubits and obtain several 
results like an improvement in the size of the stability fi- 
delity basin, analytical formulas for the distillation flows, 
phase diagrams, etc. In Sect. II Vl we apply the distilla- 
tion protocols for the purification of non-diagonal mixed 
states that arc more easily realized experimentally. In 
Sect. we study the large- Z? limit of these protocols. 
In Sect. IVll we present a detailed investigation of the re- 
lationships between quantum distillation protocols and 
renormalization methods for quantum lattice Hamiltoni- 
ans. Scction lVIIl is devoted to conclusions. In appcndixIXI 
we find the general solution for the distillation recursion 
relations used in the text in the general case of qudits. 

II. SIMPLE DISTILLATION PROTOCOLS 
WITH QUBITS 

Our starting point is the orthonormal basis of Bell 
states formed by the first qubit belonging to Alice and 
the second to Bob: 

|<I>±):=i=[|00)±|ll)] 

f (1) 
:=^[|01)±|10)] 

We shall use the word "simple" applied to the distillation 
protocols to denote that the mixed state we shall be deal- 
ing with is made up of a combination of one state in the 
set S := {|$^), |^~)} of Bell states that have coincident 
bits in Alice's and Bob's qubits, with another state in 
the set A := {|*+), |*")} of Bell states that do not have 
coincidences. Thus, wc have 4 possible combinations to 
do this type of entanglement distillation. 

To begin with, we shall choose the following mixed 
state in order to set up a simple distillation protocol 

p++:=F|<I.+)($+| + (l-F)|vI/+)(vI/+| (2) 

Alice and Bob will also need to apply the CNOT gate 
defined as usual 

UcNOTm) ■=\i)\i®j), i,J=0,l. (3) 

The distillation protocol can be arranged into a set of 
5 instructions or steps |^.[T^.p!^: 



Distillation Protocol for Qubits 

1. Set up p — > p® p with fidelity F. 

2. Apply bilateral CNOT gate: [/bcnot- 

3. Alice and Bob measure target qubits. 

4. Classical communication of results: retain coinci- 
dences {OaOb or 1a Is)- 

5. Go to step 1) with p' with fidelity F' > F. 

The simplicity of this protocol also relies on the fact 
that we do not need any depolarization step, as it is the 
case when dealing with Werner states 0. In Fig. ID we 
show a schematic picture of a single aplication of the 
purification method. Let us comment on the outcomes 
corresponding to the most relevant steps in this protocol. 
After step 1/, the 4-quit mixed state p(E)p shared by Alice 
and Bob reads as follows 

+ F(l - + |*+$+)(*+$+|] (4) 

+ (1 -F)2|1f+^'+)(^'+^'+| 

In step 2/, Alice and Bob apply bilaterally the CNOT 
gate taking their first qubit as source and their second 
qubit as target, i.e., qubits first and third are source 
qubits while qubits second and fourth are target qubits. 
To obtain the transformed mixed state we must deter- 
mine the action of the bilateral CNOT gate ?7bcnot [3 
on the states of the form \(pAfB)- The results of this 
computation are shown in Tabled With the help of this 
table we find 

C/BCNOTP++ «)p++%CNOT = -F^ | | 
+ F{1 - i^)[|$ + ^' + >($+^'+| + + + + (5) 
+ (1 + 

After steps 3/ and 4/, Alice and Bob measure their tar- 
get qubits and retain their source qubit whenever they 
find, via classical communication, the same results: ei- 
ther OaOb or Ia^b- This fact selects the state \^~^) as 
the only admissible possibility for the target state. Thus, 
only the first and third terms in the RHS of lO survive to 
this process and the resulting 2-qubit state p'_^__^_ is again 
of the same form as the original starting state Q in step 
1/, but with a higher fidelity F' > F. In fact, we get 

pV+:=F'|$+)(<i>+| + (l-F')|*+)(*+| (6) 
with the new fidelity being 

F' = (7) 

This relation defines a recursion scheme for entanglement 
purification: starting with say Np pairs of Bell states of 
fidelity F, after every application of the whole protocol 
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FIG. 1: Schematic representation of the distillation protocol 
by Alice and Bob. Originally, two pairs of shared entangled 
qubits represented by enclosed dots are transformed into a 
single pair of higher purity (doubly enclosed dots) . 
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TABLE I: This table shows the results of applying the bi- 
lateral CNOT gate to certain pairs of Bell states needed to 
distillation. 



we obtain ^ pairs of higher fidelity F' > F. Thus, pu- 
rification is achieved at the expense of halving the num- 
ber of Bell pairs. The fixed points Fc of the recursion 
relation are defined as F'{Fc) := Fc and they are 
given by Fc = 0,^,1. The fixed points Fc = 0,1 are 
stable, while Fc = i is unstable. The best way to recast 
these qualitative properties of the flow equation for the 
fidelities Q is to draw the corresponding flow diagram 
as shown in Fig. |2 



We can answer this question affirmatively by comput- 
ing the action of the bilateral CNOT gate on the tensor 
product of these mixed states®. With a similar analysis 
which has led us to Table we obtain 

UbcnotP+~ ® P+~Ubcnot = F^|$+$+)($+$+| 
+ {1- F)2|*+$-)(*+4>-| 



-I- (1 - 

C^bcnotP — ® P--C/bcnot = 
+ (1 - Ff\^+<p-)(^+^-\ 



I- |^-<If-)(<If-^-|] (9) 



We now realize that if we proceed to measure the target 
bits and classical communicate the results, we do not end 
up with the same type of mixed state as we had started 
with. That is, the protocol as it stands is not valid since 
it does not yield invariant mixed states. This problem 
has a solution provided we introduce an additional step 
prior to the measurement of the target qubits by Alice 
and Bob. This additional step corresponds to a local 
unitary operation Ua ^ Ub performed by Alice and Bob 
on their source qubits. The form of this local unitary 
depends on the mixed state we are distilling. We find 
the following results: 

Step 2'. Alice and Bob apply a local unitary transforma- 
3 Ub to their source qubits: 
Ua = lil + i)icr^ + ay), Ub = l{l-i)i<T^- ay). 



tion Ua 

For p+^, 

For p ^ 

For p 



Ub = 5(1 + cr.y). 



Ua 

Ua^o,, Ub = 1. 
After this extra step, we can guarantee that the result- 
ing 4-qubit mixed state has the appropriate Bell pairs 
at the source qubits so as to produce the same original 
state, once steps 3/ and 4/ are performed. Moreover, it 
is straightforward to prove that the new fidelity for these 
3 protocols is also given by the same recursion relation 
(0 as in the first protocol. 

Finally, if the constraint that the state of fidelity F 
must be a state is relaxed, then there are two addi- 
tional possible mixed states whose analysis can be carried 
out in a similar fashion. 



FIG. 2: Flow diagram for the fidelity F of the distillation 
protocol given by the recursion relation 0. 



III. MULTILEVEL EXTENSIONS OF 
DISTILLATION PROTOCOLS 



Next, we may wonder whether it is possible to devise 
distillation protocols for the three possible combinations 
of Bell states, namely, 

p+_ :=F|$+)(ci>+| + (l-F)|1'-)(vI/-| 

p^+:=F\^-){^-\ + {l-F)\^+){^+\ (8) 
:= F|$-)($-| + (1 - 



In order to generalize the simple distillation protocol 
of the previous section to the case of qudits, we must 
notice that the two main ingredients in that distillation 
protocol are: 

i) the CNOT gate, 

ii) the BeU states (QJ. 

Regarding the CNOT gate, the extension of this gate 
to deal with qudits is not unique. As has been noted in 
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[T3,[ll, the CNOT gate for qubits ® has 3 properties 
that make it special, namely 



U. 



CNOT 

t 



^CNOT' 
^CNOT = C^CNOT, 

i® j =^0 ^ i = j. 



(10) 



The extension of the CNOT gate for audits that satisfies 
these 3 properties (fTn|l is given by 



f/cNoxN)!.?) N)Nej), = 0,. 



(11) 



where iQj := i—j, mod D, denotes substraction modulus 
D. This is the definition that we shall adopt throughout 
this paper, unless otherwise stated. 

As for the higher-dimensional extension of Bell states 
O, we shall also take the following generalization 

\^k,) :=C/cNOT[(C/F|fc))® li)], fc,j=0,...,i^-l (12) 
where C/p is the quantum Fourier transform (QFT) 



UT?\k) := 



1 



D-l 



'■\y), 



(13) 



which reduces to the Hadamard gate when dealing with 
qubits {D = 2). As a matter of fact, we can readily check 
that for the special case of qubits D = 2 we recover the 
standard Bell pairs in the following form 



l^-oo) = 1$' 



),l*Ol) = 1*^ 



(14) 



Moreover, using the generalized CNOT gate Hll|) . the 
generalized Bell states are given by 



D-l 

E 

y=o 



c " \y)\yej) 



(15) 



With these extensions of the CNOT gate and the Bell 
states, we can set up a generalization of the simple dis- 
tillation protocols of Sect. II for qudit states. These 
protocols have the same 5 steps as before. 
Step 1 /. We shall assume a general diagonal mixed state 
of the form 



Z5-1 
fcj=0 

1 =: E ^'^J' 

k,j=0 



(16) 



where qkj are normalized probabilities. For non-diagonal 
mixed states, we refer to Sect. IV. Then, Alice and Bob 
share pairs p x p oi these states (|16ll . 
Step 2/. Alice and Bob apply bilaterally the generalized 
CNOT gate (|llf) . To know the result of this operation 



on the state (fTB|) we need a previous result about the 
action of the gate ?7bcnot on pairs of generalized Bell 
states p2|l . After some algebra, we arrive at the following 
expression 

C/BCNOT|*fcj)|«'fe'j') = |«'feffife'j)|^'Defc'je/> (17) 

This is a fundamental result for it means that the space 
of two-pairs of generalized Bell states is invariant under 
the aetion of the generalized bilateral CNOT gate. This 
is a very nice result that condenses in a single formula 
all the possibilities for the outcome of the action of the 
CNOT gates on Bell states, in particular, the whole table 
employed by Bennett et al. in |3| for the case of qubits is 
contained in equation 117|) . This property is essential in 
order to have a closed distillation protocol. Actually, it 
would have been enough to have obtained only the source 
qubits as generalized Bell states. 

Then, with the help of this property (|17|l we obtain 
the action of C/bcnot on pairs of p states, as follows 



C^BCNOT/O <Xi pC^BCNOT = ^ ^ qkek',]qk'j' 
kj=0 k',j'=0 

\^kj-^Dek',jer){'^kj'ifDek',jej'\ 



(18) 



We see that this state is already of the same form in the 
source qubits as the original p 

Step 3/. Alice and Bob measure their target qubits in 
(|18|l . To see the result of this measurement, let us write 
the explicit form of the target qubits, namely 



1 



"^Dek'^jer = 7^ E e-^|z)|ze (je/)). (19) 

Therefore, coincidences between Alice's and Bob's target 
qubits will happen only when the following condition is 
satisfied 



z = zQ{jQj')-^j=j'. 



(20) 



Step 4/- After their measurement, Alice and Bob com- 
municate classically their result so that they retain the 
resulting source Bell pairs only when they have coinci- 
dences, and discard them otherwise. The resulting net 
effect of this process is to produce a Kronecker delta 
function Sjji in the target qubits. More precisely, the 
resulting unnormalized mixed state is given by 

D-l D-l 

P' - E 9fcefc'j%'/'5jy|*fcj)(*fc,|. (21) 

fc,i=0 k',j'=0 

Therefore, we end up with a diagonal mixed state of the 
same form as the starting one 



D-l 

^ q',,\^k,){^k,\ 
k,j=0 



(22) 
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with the new probabihties given by 

2^fc,7=o /^k'=Q Qkek'jQk'] 



(23) 



This is a generahzed recursion relation that inchides eq. 
Q as a particular instance. 

Step 5/. Alice and Bob start all over again the same 
process with the initial state now being p' in (|22I23() . 

The nice feature of these generalized distillation pro- 
tocols for dealing with qudits is the fact that we have 
at our disposal explicit analytical formulas H23|l for the 
evolution (flow) of the different weights (probabilities) 
of the generalized mixed states to be purified. As these 
distillation protocols are too general, it is worthwhile to 
consider some particular cases of interest separately. The 
general solution to the distillation recursion relations H23|) 
is presented in the appendix \K\ Wc hereby provide the 
following analysis of some examples: 

i). Let us investigate the closest generalization of the 
simple protocols introduced in Sect. II. Thus, let us con- 
sider the following type of initial mixed state 



M 

p■.= Y,Q^\'^0^){'^0^\, M<D-1, 
1=0 
M 



(24) 



This corresponds to working with the subset of all possi- 
ble generalized Bell states of the form {|^E'oi)}iLo^- Inter- 
estingly enough, this includes the case of the state p++ 
in for = 2. The recursion relations (|23|1 for this 
special subset of states takes the following simpler form 



/ 

9t 



if 



(25) 



For M — 2, i.e., considering a mixed state formed of 
just two Bell states of the form |^oi)j the protocol has 
the following recursion relation 



qf + (1 - * 



(26) 



where here the index i stands for any possible pair of Bell 
states of the type l^'oi)- In other words, we have found 
a direct Z?-dimensional generalization of the distillation 
protocols for qubits in Sect. II, with qi := F. 

a). For M = D — \ and taking go F and qi :~ 
i = 1, . . . , Z) — 1 we can find a more advantageous 
protocol than the previous one. In fact, in this case we 
find that 



D-l 



F2 



F2 



Z5-1 



(27) 



The fixed points of this recursion relation are now given 
by Fc = 0, -^,1. Despite being a non- linear recursion 



F' 




FIG. 3: Plots of the distilled fidelities F' as a function of the 
original fidelity F for several values of the dimension D of the 
qudits: D = 2, 4, 10. 



relation, 127|l admits an explicit analytical solution for 
the general term of the series Fk given by 



F,: - 



F^-^") + {D - I) 



l-F 



D-1 



k>l,Fo:^F. (28) 



From this solution, wc inmediately find that the fixed 
points Fc = 0, 1 are stable while Fc = is unstable. 

In Fig. 13 we plot the function F' = F'{F) for several 
values of the dimension D. From the analysis of these 
curves we inmediately obtain the corresponding flow di- 
agram that we represent in Fig. 21 We check that for 
D ~ 2 we recover the flow diagram corresponding to 
standard qubits (Fig. ^ . 



D 



Fc 



FIG. 4: Flow diagram for the fidelity F of the generalized 
distillation protocol for qudits given by the recursion relation 

We sec from Fig. 0| that the stability basin is increased 
with respect to the case of standard qubits, as in Fig. |21 
This means that we can start with a mixed state having 
a fidelity F with respect to the Bell state l^'oo) lower 
that i and we still will succeed in purifiying that state 
towards fidelity close to 1. Thus, we have found that it 
is more advantageous to distill a given Bell state 1^*00 ) if 
we prepare the mixed state p in H24|) in the form 



D-l 



p := F|1'oo)(1'oo| + jy—^ l*0')(«'o^|, 



rather than using just one single of those states 
p := F|*oo)(«'oo| + (1 - J^)|*o«)(*Oz|, i ^ 0. 



(29) 



(30) 
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We may wonder how is it likely for Alice and Bob to 
obtain the same values (coincidences) after measuring the 
target qudits in the step 3/ of the distillation protocol. 
Let us denote by Vab this probability which will depend 
on the value F of the fidelity. From equations (|21|l and 
()27|l we find this probability of coincidences to be 



Vab{F) = F'' 



D - 1 



(31) 



The minimum of this probability is at i^o = ^i^d its 
value is Vab{-j^) = j)- Likewise, 7^ab(1) — 1- Thus, we 
find that the probability is lower and upper bounded as 

One is also interested in knowing the number of steps 
iC(e, Fq) needed to achieve a certain final fidelity close 
to 1, say 1 — e, starting from an appropriate initial fi- 
delity Fq > j^. This number can be computed from our 



analytical solution H28|) from the condition 



(32) 



Thus, we find the following analytical formula for the 
number of steps needed to obtain a certain degree of fi- 
delity e as a function of the initial fidelity _fo > "^ji-e.. 



K{e,Fo) 



l0g2 



{l -e){D- l) 
1-F„ 



it) 



(33) 



In Fig. |31 we plot the number of iterations (|33|l for a 



K 




FIG. 5: Plot of the number of iterations K{e, Fa) lESJ to 
achieve a final fidelity of F = 0.99 as a function of the initial 
fidelity Fq and for several values of the dimension D of the 
qudits: D = 2, 4, 10. 

given value of the final fidelity 1 — e that we take as the 
fixed value of 0.99, and then we find how is the depen- 
dence on the initial fidelity Fq. We see that for a given 
admissible value of Fq, the lowest number of iterations 
corresponds to the protocol with the higher value of the 
qudits dimension D: 




FIG. 6; The desired fidelity q'q for amplification after a single 
application of the distillation protocol for qutrits 1361 as a 
function of the previous fidelitites qo,qi. 



in). For qutrits, D = 'i, the most general diagonal 
mixed state with the allowed Bell states taking values on 
the set {|*oi)} is 



p := (7o|*oo) +f?i|*oi) 
1 =: qo + qi + q2 



'721* 



02/ 



(35) 



Let us assume that the state we want to purify is l^'oo)- 
Now, our recursion relation for our fidelity go depends on 
two variables, namely, 



Qo = 



qI 



qI 



2 ' (i-go-gi)2' 



(36) 



K{eo,Fo)D^ > K{eo,Fo)D„ for Di > D^. 



(34) 



and a similar equation for qi with q^ qi. In Fig. |H| 
the dependence of the function fidelity q^ = qQ{qo, qi) for 
qutrits is plotted. We observe that it is a monotonous 
incresing function which gurantees that the initial fidelity 
will flow towards 1, under certain conditions. To find 
these conditions, wc find that the set of fixed points of 
these recursion relations is given by 

(go,9i)c = {(0,0),(i,0),(0,i),(i,i),(i,i),(l,0),(0,l)}. 

We have also found the flow diagram associated to these 
recursion relations which is now two-dimensional and 
we show it in Fig. [3 From this diagram we see that 
the purification protocol is successful in arriving to the 
maximum fidelity qa ~ 1 provided the initial fidelity 
lies in the stability basin of the fixed point (1,0) which 
is given by the trapezoid formed by the set of points 
(i,i),(i,i),(i,0),(l,0). 



IV. DISTILLATION OF NON-DIAGONAL 
MIXED STATES 

So far, we have been investigating the properties of 
distillation protocols applied to mixed bipartite state of 
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(0,1) 




(1,0) 



FIG. 7: Two-dimensional flow diagram associated to the dis- 
tillation protocol for qutrits 1361 . 



diagonal form such as those in H22|) . (|23|) . It is apparent 
that once we have a general result for the operation of 
the C/bcnot gate of generalized Bell states (|17ll . we can 
also deal with non-diagonal mixed states, namely, 



D-l D-l 

k,j=0 fc' j'=0 
D-l D-l 

1 =: X! Ikjk'f, qkjk'j' > 0. 

k,j=0 k',j'=0 



(37) 



Since this is a too much general state, we prefer to extract 
from this class of non-diagonal states one type which we 
believe it may have potential applications. 

Let us imagine that Alice and Bob are manipulating 
bipartite qudit states that are diagonal in the computa- 
tional basis. More explicitly, the entangled state they 
want to purify is of the form 



l*a):^^El«> 



(38) 



while states which are non-diagonal are considered as 
acting as disturbing noise that they want to get rid of. 
Specifically, this noise will be represented by the state 



D-l 



(39) 



Then, in order to achieve their goal of purifying states 
of the diagonal form \'^d) with respect to non-diagonal 
states l^'o), they set up a distillation protocol based on 
sharing copies of the following mixed state 



J^|*d)(*d| + (l-i^)|4'o)(*o|. 



(40) 



We envisage that this scenario is physically feasible since 
we can imagine that the computational basis is realized 



in terms of some physical property taking values on z = 
0, 1 and that Alice and Bob have a mechanism to 

select when they have \ii) coincident qudits (or diagonal) 
from \ij) non-coincident qubits (non-diagonal). 

To proceed with the distillation of the state p in H40|l . 
we first must express the states 1"^^) and |5'o) in the basis 
of the generalized Bell states, with the result 



|*d) = 1*00 

I*. 



D-l 



(41) 



Next, Alice and Bob share two pairs of non-diagonal 
mixed states 

/9(g) p = F2|1'oo«'oo)(*oo«'oo| 
F(\ — F) ^ l 

+ J2 [l*00*0^>(*00*Oil + |*0j*00)(*Qj*00|] 



(l-F) 



D-l 



(D-l) 



J2 l*0..*0fc)(«'0j*0/|, 



(42) 



and they apply bilaterally the CNOT gate to it H17|l with 
the result 

%cnotP «) P%CNOT = -F'^|*oo*oo)(*oo*ool 
Fil-F) 



D 



[|*oo«'o,e«)(*oo*o,ej| + l^-o.^oz) (^'oi^ojl] 



i.j=i 



' ' ij,k,l = l 



(43) 



The process of measuring the target qudits and retaining 
the source qudits when upon classical communication Al- 
ice and Bob find coincidences in their measures amounts 
to retaining the terms in H43|) that have the state IvPoo) 
in the target qudits. This means that only the first term 
and part of the last term in H43|) contribute to the fi- 
nal source mixed state, which takes the following form 
without normalization 



f2|«'oo)(*oo| + ^^-^|*o)(*o| 



(44) 



Upon normalization, we arrive again at a non-diagonal 
mixed state of same form as the one we started with 
p' = F'|*d)(*d| + (1 - i^')l*o)(*o|, but with a new 
fidelity F' given by 



F' = 



F2 



D-l 



(45) 
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Let us notice that this is precisely the same recursion 
relation that we found in Sect. Ill in a different context 



(X) 



CONTINUUM LIMIT OF QUDIT 
PROTOCOLS 



For the general case represented by the recursion re- 
lations H25() we can also find the general solution for the 

fc-th iteration ql''\ i = 0, . . . , D — 1 starting from their 

initial values qf^ satisfying J^Z^^ 1- ^^^^^ 

following solution 




(fc) 
9i = 



.(0) 



E 



M 



(0) 



(46) 



Let us assume that the maximum initial value is M := 
max {(/I"''} and it is p times degenerate. Then, using 
the general solution (|46|l we can inmeditely find the fixed 
points after the evolution with the recursion relations. 
We find 



lim q, 

k — ^oo 



(fe) 



0, if < M, 



lim = - , if q'"' = M. 

k^oo p 



(47) 



From the analysis of these distillation protocols and the 
way they operate we arrive at the conclusion that they 
resemble a sort of amplitude amplification quite similar 
to what happens in the Grover algorithm where there 
exists what is called quantum amplitude amplification. 
However, there is an important distintion between both 
procedures: in the distillation method, the maximum am- 
plification is attained asymptotically, while in Grover al- 
gorithm it is achieved periodically. The reason for this 
difference relies on the fact that the distillation process 
is not unitary (since we make measurements and discard 
states), while Grover is unitary. Thus, we propose to 
refer to the distillation protocol as weight amplification, 
since it is certain probability weights of the intial mixed 
states, and not amplitudes, what are being amplified. 

When D is very large, we can approximate the prob- 
ability weights qf^"^ taking values on the discrete set 
{0, 1, . . . , Z? — 1} 3 i, by a density function q{xY^'> de- 
fined on the real interval [0, 1]. This is achieved by intro- 
ducing the variable x G [0, 1] defined as x := lAx with 
Aa; jJ , . Thus, in the limit D oo{Ax 0), we get 



a probability density as ql^\iAx) q^'^\x)dx. It is also 
normalized as 



q^°\x)dx = 1. 



(48) 



Likewise, we can take the continuum limit of the general 
recursion equation H46|) in order to obtain the probability 



FIG. 8: Evolution of the probability density under the itera- 
tive application of the distillation protocol in the continuum 
limit 14911 . Starting with a parabolic distribution fc = 0, we 
show the resulting profiles after steps fc = 1, 2, . . . , 6. 



density q{x)^^'^ after k steps of the distillation protocol. 
This is given by 



gW(a;) 



(49) 



This is a closed analytical equation that provides us with 
the evolution of the probability density for any initial 
probability profile q^'^\x). In Fig.|Slwe plot this evolution 
for an initial distribution of a parabolic form q^*^^ [x) = 
6(.T — x"^). We see how as we increase the step k of the 
distillation, the new distributions get peaked around the 
highest value of the initial distribution, which is a; = ^ in 
this particular case. This behaviour illustrates the idea 
of the weight amplification and is in agreement with the 
results ()47|l for the fixed points of the flow equations. 



VI. QUANTUM DISTILLATION AND 
QUANTUM RENORMALIZATION 

It is interesting to notice the analogy between the re- 
cursive distillation process represented by the equation 
O and Fig. [3 and the truncation process in the Renor- 
malization Group analysis of certain quantum lattice 
Hamiltonian models, specifically, the ITF model (Ising 
in a Transverse Field) 131, H^- The basic idea of a QRG 
method is: i/ elimination of high energy states plus, ii/ 
iterative process. This is precisely what happens in a 
quantum distillation process which wc have seen in the 
precceding sections, achieving a purification of a mixed 
state by means of discarding states and a recursive pro- 
cedure. This relationship can be made even closer if 
we briefly recall what a quantum renormalization group 
(QRG) method is. The subject of the distillation is a 
mixed state operator p, while that of the renormaliza- 
tion is a quantum Hamiltonian operator H . A summary 
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C^Uciritum J_/istill&tioii 


LldllL Lllil XVvjr 


Mixed State p 


Quantum Hamiltonian H 


Computational Basis 


Local Site Basis 


Bell Basis 


Energy Basis 


Alice & Bob Tensor Product 


Blocking Method 


L.O.C.C. 


Truncation Operator 


Maximum Fidelity 


Minimum Energy 


RG-Flow Diagram 


Distillation-Flow Diagram 



TABLE IL Summary of the comparative analysis between the quantum distillation process and a quantum renormalization 
group method for lattice Hamiltonians. 



Initial chain with 8 sites 



FIG. 9: Block decomposition of the Heisenberg chain in 3-site 
blocks. 



of these relations is presented in Table ^1 that will be 
deduced along the way. The easiest way to present 

Fo Fo Fo Fo Fo Fa Fo Fo 

;h: s s s s s s s 

X Fi Fi y X -Fi ^ X Fi ^ 

8 ft n n 



J.-2 

n 



Fa 



FIG. 10: Example of distillation process starting with 8 pairs 
of mixed states for Alice (and the same amount for Bob). 
After 3 steps the original fidelity Fq is improved up to a final 
value of F3 (we assume full success for simplicity). 



the QRG method is with an example of quantum lattice 
Hamiltonian like the isotropic Heisenberg model on a ID 
chain: 



N 



H = J^S, • S 



QRG 




FIG. 11: Example of renormalization QRG process for an 
initial chain with 8 sites in 3 steps (making blocks of 2 sites 
each) . 



of the distillation process, as shown in Fig. ^| This is 
to be compared with the similar iterative process in the 
QRG method in Fig. CT] 

The block Hamiltonian is then 



i^B = J(Si • S2 + S2 • S3) 



82 + 83)2 



- (81 



83)^1 



(51) 



The label B here stands for Block and not for Bob. 
The diagonalization of Hb is straightforward using the 
Clebsch-Gordan decomposition of the tensor product of 
3 irreducible representations of spin 5 = ^, 



(52) 



(50) In particular, the ground state (GS) is given by 



with Si spin-i operators at site i of the chain. The local 
site basis {| t)} corresponds to the computational 
basis {|0), |1)}. Much like this latter basis is not enough 
for doing the distillation, the local site basis needs to be 
complemented with another type of basis. To see this, 
let us start the KG process with the block decomposition 
of the chain in blocks of ub = 3 sites as shown in Fig. 
This blocking method in QRG corresponds to the tensor 
product of Alice and Bob's shared states at the begining 



|^)Gs = ^[2|TiT)-UTT)-|TTi)], (53) 

which is a spin doublet (with a similar exprcsion for the 
other state |-1J-)gSj with the spins reversed). This fact is 
peculiar of the 3-site block and it is the main underlying 
reason for using a block of that size in the QRG (this fact 
is model dependent: for the ITF model, the blocking is 
with Ub ~ 2 sites |23,|23, Fig. In the energy basis, 
the block Hamiltonian is diagonal and this corresponds 



10 




FIG. 12: QRG renormalization of Hb- 



to the Bell basis for the mixed state p in the distillation 
process. 

Now, the truncation of states amounts to retaining the 
state of lowest energy (doublet) and discarding the re- 
maining 2 excited states. This reduction scheme is of 
the form 2'^ = 8 — > 2. This truncation corresponds 
to discarding unwanted states of non-coincidences in the 
distillation process. The new effective site is again a spin- 
^ site as shown in Fig. 1121 The RG-truncation is imple- 
mented by means of a truncation operator O constructed 
from the lowest energy eigenvalues of Hb retained during 
the renormalization process. In this example, O is con- 
structed from the lowest energy doublet in the Clcbsch- 
Gordan decomposition (|52|) . namely H53|l . Similarly, in 
the distillation process wc have that the tensor product 
of Alice and Bob's states can be decomposed into states 
with coincident qubits in the target, denoted hy (0J, 
and states with non-coincident qubits in the target, de- 
noted by Q' i-^-' 

pAB®PAB^{Y.Pc)®C^P^c)^ (54) 

a a 

where the sum in a runs over a certain number of mixed 
states of 4 parties. Notice that this stage is similar to the 
RG-stage represented by equation ifS^ . Next, an elimi- 
nation process similar to the RG-truncation is performed 
by means of LOCC operations (measurements and classi- 
cal communication) that retains only the bipartite states 
embedded in the Pq states. 

Then, the renormalization of the block Hamiltonian is 
simply 

Hb' = OHbO^ =Eo = - J. (55) 

Similarly, we could have expressed the state-elimination 
of the distillation in previous sections in terms of an trun- 
cation operator, say Od, such that the new mixed state 
p'j^B is obtained as 

Pab = Od{pab ® pab)0'^d- (56) 

In the case of the quantum Hamiltonian, we still need 
extra work since there are interaction links between 
blocks (see Fig. O. These are absent in the distillation 
protocol. However, the renormalization of the interblock 
Hamiltonian Hbb follows also the same prescription as 
in H55|) and we arrive at 

RG 2 

JS" • S; * '^^3-' ■ ^n-f-l' (57) 



1/ Block Decomposition: H = Hb + Hbb- 

2/ Diagonalization of Hb- 

3/ Truncation within each Block: O- 

4/ Renormalization: Hg = OHbO'' , H'bb = OHbbO 

5/ Iteration: Go to 1/ with H' ^ H'g + H'sg. 



TABLE III: Steps of the quantum renormalization group 
method (QRG) for lattice Hamiltonians. 



n n + 1 




n n + 1 



FIG. 13: QRG renormalization of the interblock Hamiltonian 
Hbb- 



where wc denote by n and n -\- \ two successive blocks 
in the original lattice (see Fig. 113(1 that become two suc- 
cessive sites (see Fig. I12f) in the new lattice after the 
renormalization. We can collect all these steps in Ta- 
ble |llTl [111 , [H. This table should be contrasted with 
the similar table for the distillation process that can be 
formed with the steps explained in Sect. ITU 

The outcome of the RG-method is that we obtain the 
correct RG-flow for the coupling constant J — > 0, sig- 
nalling a gapless system plus an approximate estima- 
tion for the ground state energy, which by means of the 
variational principle, it is an upper bound for the ex- 
act energy. Therefore, the QRG is an energy minimiza- 
tion procedure. Likewise, the purification process pro- 
duces a protocol for fidelity maximization along with a 
distillation-flow diagram. This completes the relation- 
ship established in Table ^] between quantum distillation 
and quantum renormalization. 

VII. CONCLUSIONS 

The field of quantum distillation protocols has become 
very active in the theory of quantum information due to 
the central role played by entanglement in the quantum 
communication procedures and its tendency to degrada- 
tion. 

In this work we have been interested in several exten- 
sions of the purification protocols when dealing with mul- 
tilevel systems (qudits) instead of the more usual qubit 
protocols. We have seen the various advantages of having 
distillation methods for qudits systems as compared with 
the simple case of qubits. We have also obtained the gen- 
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eral form of the solution to the distillation recursion rela- 
tions and several particular solutions have been studied 
explicitly. We have developed the relationship between 
quantum distillation protocols and quantum renormal- 
ization group methods, something which is interesting in 
itself and could serve as a guide for possible extension of 
purification methods. 

We would like to mention that the possibility of work- 
ing with qudits systems has become quite realistic in the 
recent years. For instance, it is possible to realize mul- 
tilevel systems in terms of the orbital angular momen- 
tum of photons, instead of the more standard polariza- 
tion (qubit) degree of freedom [l^, H, [2i|. Yet another 
possibility is to use the so called multiport beam splitters 

[23, ia, lal, 0. 

There are several ways in which this work can be ex- 
tended. One is the consideration of noise as as source of 
errors during the distillation protocol itself. Another one 
is to allow the possibility of having these distillation pro- 
tocols for qudits be embedded into a quantum repeater 

protocol Ea, [3. 

Acknowledgments. This work is partially supported by 
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Thus, are unnormalized probability weights. The re- 
cursion relations they satisfy can be read as follows ljAl|l : 
for a fixed second index j, the unnormalized weights g^""* 
at the step n of the distillation process are obtained as the 
convolution over the first indices k of the unnormalized 
weights g^" in an earlier step. This fact calls for the 
introduction of the Fourier transform in order to analyze 
the relations (|A1|) . Let us introduce the new variables 
defined as 

^^-E^'^^i?. (A3) 

Now, using the properties of the convolution and the 
Fourier transform it is inmcdiatc to arrive at a simpler 
recursion relation 

= Ti?!"-!)] ' , (A4) 
which can be iterated all the way down to the initial step 



APPENDIX A: GENERAL SOLUTION OF THE 
DISTILLATION RECURSION RELATIONS 

In this appendix, we look for more general solutions 
to the general distillation recursion relations (|23|) than 
those studied in section ITTll To this end, it is convenient 
to introduce auxiliary variables ^["^ defined by 

D-l 

Jn) _ V- (n-1) (0) _ (0) , ... 

^kj ^ 9kek'j9k'j ' Skj ~ Ikj ' y^'^) 



k'=0 



SO that the real weights g^"' arc related to these auxiliary 
variables as 



(n) 
Ik/ 



(") 

9 k/ 



l^Li=0 yii 



k] 



. kj 



(2") 



(A5) 



Fourier transforming back to the unnormalized variables, 
we get 



(") 
9 k/ 



1 

D ^ 



fc=0 



D-\ 



- (0) 

"Ik/ 



(2") 



(A6) 



(A2) 0. 



from which we also obtain the normalized probability 
weights q//^ upon normalization (|A2p . In particular, for 
the case of qubits treated in Scct.^ D — 2 and if we also 
restrict ourselves to weights of the form qqq = q^i = 
1 — F, qio = qii = 0, we again obtain from the general 
solution ljA6p the simple recursion relation in equation 



[1] M. Lewenstein, D. Bruss, J.I. Cirac, B. Kraus, M. Kus, J. 
Samsonowicz, A. Sanpera, R. Tarrach, "Separability and 
distillability in composite quantum systems -a primer-" , 
J. Mod. Opt. 47, 2841 (2000). 

[2] D. Bruss, J.I. Cirac, P. Horodecki, F. Hulpke, B. Kraus, 
M. Lewenstein, A. Sanpera, "Reflections upon separabil- 
ity and distillability" , J. Mod. Opt. 49, 1399-1418 (2002). 

[3] C.H. Bennett, G. Brassard, S. Popescu, B. Schumacher, 
J. A. Smolin, W.K. Wootters, "Purification of Noisy En- 
tanglement and Faithful Teleportation via Noisy Chan- 
nels", Phys. Rev. Lett. 76, 722-725 (1996). 

[4] J.W. Pan, C. Simon, C. Brukner, A. Zeifinger, "Feasi- 
ble Entanglement Purification for Quantum Communi- 
cation", Nature 410, 1067 (2001). 



[5] N. Gisin, "Hidden quantum nonlocality revealed by local 
filters", Phys. Lett. A210, 151, (1996). 

[6] P.G. Kwiat, S. Barraza-Lopez, A. Stefanov, N. Gisin, 
"Experimental entanglement distillation and 'hidden' 
non-locality," Nature 409, 1014 (2001). 

[7] P.W. Shor, "Scheme for reducing decoherence in quan- 
tum computer memory", Phys. Rev. A52, 2493-2496, 
(1995). 

[8] A.M. Steane, "Error correcting codes in quantum the- 
ory", Phys. Rev. Lett. 77, 793, (1996). 

[9] D. Gottesman, "A theory of fault-tolerant quantum com- 
putation," Phys. Rev. A 57, 127-137 (1998). 
[10] J. Preskill, "Fault-tolerant quantum computation." 
.quant-ph/9712048. 



12 



C. Macchiavello, 
QPA procedure" . 

C. Macchiavello, 

D. Bouwmeester, 



C. H. Bennett, D.P. DiVincenzo, J. Smolin, W.K. Woot- [23 
ers, "Mixed-state entanglement and quantum error cor- 
rection", Phys. Rev. A54, 3824 (1996). 

D. Deutsch, A. Ekert, R. Jozsa, C. Macchiavello, S. [24 
Popescu, A. Sanpera, "Quantum Privacy Amplification 

and the Security of Quantum Cryptography over Noisy 
Channels", Phys. Rev. Lett. 77, 2818-2821 (1996). 

"On the analytical convergence of the [25 
Phys. Lett. A246, 385 (1998). 
"Quantum Privacy Amplification", in 
A. Ekert, A. Zeilinger (eds.). The [26 
physics of quantum information, Springer- Verlag 2000. 
W. Diir, H.-J. Briegel, J.I. Cirac, P. ZoUer, "Quantum re- [27 
peaters based on entanglement purification", Phys. Rev. 
A. 59, 169 (1999). 
W. Diir, H.-J. Briegel, J.I. Cirac, P. ZoUer, "Quantum [28 
repeaters for communication", quant- ph/9803056 
M. Horodecki, P. Horodecki, "Reduction criterion of sep- 
arability and limits for a class of distillation protocols", 
Phys. Rev. A59, 4206-4216, (1999). [29 
G. Alber, A. Delgado, N. Gisin, I. Jex, "Generalized 
quantum XOR-gate for quantum teleportation and state 
purificationin arbitrary dimensional Hilbert spaces", [30 
quant-ph/000802. 

G. Alber, A. Delgado, N. Gisin, I. Jex, "Efficient bipar- 
tite quantum st ate purification in a rbitrary dimensional [31 
Hilbert spaces" , |quant-ph/0102035| 

H. Bechmann-Pasquinucci, A. Peres, "Quantum Cryp- 
tography with 3-State Systems", Phys. Rev. Lett. 85, [32 
3313-3316, (2000). 

H. Bechmann-Pasquinucci, W. Tittel, "Quantum cryp- 
tography using larger alphabets", Phys. Rev. A 61, 
062308 (2000). [33 
M. Bourennane, A. Karlsson, G. Bjork, "Quantum key 
distribution using multilevel encoding" , Phys. Rev. A 64, 
012306 (2001). 



C.H. Bennett, H.J. Berstein, S. Popescu, B. Schumacher, 
"Concentrating partial entanglement by local opera- 
tions", Phys. Rev. A53, 2046 (1996). 
J. Gonzalez, M.A. Martm-Delgado, G. Sierra, A.H. Voz- 
mediano. Quantum Electron Liquids and High-Tc Super- 
conductivity, Lecture Notes in Physics, Monographs vol. 
38. Springer- Verlag 1995. 

M.A. Martm-Delgado, G. Sierra; "Analytic Formulations 
of the Density Matrix Renormalization Group"; Int. J. 
Mod. Phys. All, 3147-3174 (1996). 

S.R. White, "Density-matrix algorithms for quantum 
renormalization groups" Phys. Rev. B48, 10345 (1993). 
A. Mair, A. Vaziri, G. Weihs, A. Zeilinger, "Entangle- 
ment of Orbital Angular Momentum States of Photons" , 
Nature 412, 3123-316 (2001). 

H.H. Arnaut, G.A. Barbosa, "Orbital and Intrinsic An- 
gular Momentum of Single Photons and Entangled Pairs 
of Photons Generated by Parametric Down-Conversion" 
Phys. Rev. Lett. 85, 286 (2000). 

S. Franke- Arnold, S.M. Barnett, "Two-photon entangle- 
ment of orbital angular momentum states", Phys. Rev. 
A65, 033823 (2002). 

M. Reck, A. Zeilinger, H.J. Bernstein, P. Bertani, "Ex- 
perimental realization of any discrete unitary operator" , 
Phys. Rev. Lett. 73, 58-61 (1994). 

M. Zukowski, A. Zeilinger, M.A. Horne, "Realizable 
higher-dimensional two-particle entanglements via multi- 
port beam splitters", Phys. Rev. A 55, 2564-2579 (1997). 
G. Molina- Terriza, J. P. Torres, L. Torner, "Management 
of the Angular Momentum of Light: Preparation of Pho- 
tons in Multidimensional Vector States of Angular Mo- 
mentum", Phys. Rev. Lett. 88, 013601 (2002). 
J. Leach, M.J. Padgett, S.M. Barnett, S. Franke- Arnold, 
J. Courtial, "Measuring the Orbital Angular Momentum 
of a Single Photon", Phys. Rev. Lett. 88, 257901 (2002). 



